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ASYMPTOTIC ANALYSIS OF FORWARD PERFORMANCE PROCESSES IN 
INCOMPLETE MARKETS AND THEIR ILL-POSED HJB EQUATIONS 


MYKHAYLO SHKOLNIKOV, RONNIE SIRCAR, AND THALEIA ZARIPHOPOULOU 

Abstract. We consider the problem of optimal portfolio selection under forward investment per¬ 
formance criteria in an incomplete market. The dynamics of the prices of the traded assets depend 
on a pair of stochastic factors, namely, a slow factor (e.g. a macroeconomic indicator) and a fast fac¬ 
tor (e.g. stochastic volatility). We analyze the associated forward performance SPDE and provide 
explicit formulae for the leading order and first order correction terms for the forward investment 
process and the optimal feedback portfolios. They both depend on the investor’s initial preferences 
and the dynamically changing investment opportunities. The leading order terms resemble their 
time-monotone counterparts, but with the appropriate stochastic time changes resulting from av¬ 
eraging phenomena. The first-order terms compile the reaction of the investor to both the changes 
in the market input and his recent performance. Our analysis is based on an expansion of the un¬ 
derlying ill-posed HJB equation, and it is justified by means of an appropriate remainder estimate. 


1. Introduction 

This paper analyzes the optimal portfolio selection problem under forward investment criteria in 
incomplete markets. Incompleteness stems from the presence of imperfectly correlated stochastic 
factors that affect the dynamics of the traded assets. Such factors have been widely used in the 
literature and model an array of market inputs, like, among others, stochastic volatility, stochastic 
interest rates, predictability of asset returns, and various macroeconomic indicators. Herein, we 
consider a pair of such factors, which are taken, however, to move at different time scales. 

The mathematical formulation of the problem of optimal investment in continuous time was 
pioneered by Merton in |Mel| . [Me2| and is usually referred to as the Merton problem. In the 
classical Merton problem the investor faces a complete market and seeks an investment portfolio 
that optimizes her expected utility from wealth acquired in the investment process. Therein, the 
investor’s utility function (or, equivalently, her preferences) is determined ex ante and does not 
change over time. The Merton problem has been studied in a variety of frameworks and we refer to 
the books |Du| . |KS| for excellent accounts of the classical results. However, the setup of the Merton 
problem has two inherent drawbacks: 1) the investor has to decide about her utility function ex 
ante and cannot adapt it to market observations; 2) the investments over different time horizons are 
typically inconsistent with each other: for example, for 0 < Ti < T 2 the solution to the investment 
problem for the time period [0,Ti] is in general not the restriction to [0, Ti] of the solution to 
the investment problem for the time period [0, T 2 ] with the Ti utility function also imposed at 
T 2 . In particular, it is not clear how to invest in a time-consistent fashion over rolling horizons 
(deterministic or random). 

Forward investment performance criteria were introduced and developed in [MZlj and [MZ2j . 
and provide a complementary setting to the traditional expected utility framework. They allow 
for dynamic adaptation of the investor’s risk preferences given how the market conditions change 
and, also, take into account the updated performance of the implemented strategies. The forward 
performance process U{t,-), t > 0 is a stochastic process adapted to the filtration of the investor 
with the properties that with probability one all functions x U{t, x) are increasing and concave 
(and, thus, can serve as utility functions); for every self-financing strategy vr and the corresponding 
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portfolio value process X'^ the process U{t,X'^{t)), t > 0 is a super martingale in the filtration of 
the investor; and there exists a self-hnancing strategy vr* such that the process U{t,X'^ (t)), t>0 
is a martingale in the filtration of the investor. The pair encodes how the preferences of 

the investor and her optimal investment decisions jointly evolve in time from the given C/(0, •). 
For a more detailed description of these criteria, further motivation and construction of concrete 
examples, we refer the reader, among others, to [E1M| . |KOZ| . |MZ5] . |MZ6] . [NT] , and [NZl] , 

The specification of the initial datum 17(0, •) is a central issue and a topic of current research in 
the forward portfolio selection approach. From the theoretical point of view, the main question is 
the characterization of the set of admissible initial data that give a well-defined solution. This was 
addressed in |MZ4j for the case of time-monotone forward processes and we refer to these results 
repeatedly herein. 

From the practical point of view, the question is how to use and translate investor’s targets 
- for example, desired upcoming investment performance or personal market views, etc. - to a 
mathematical input. Questions of this nature were studied in [MZ4| (section 5.2) where it is 
shown how an investment target on the upcoming average return can be used to infer the initial 
utility input. For example, one can think of a client presenting a fund manager with the desired 
investment target (e.g. 5% above the S&P 500 performance) and a band around the investment 
target (e.g. 4-6% above the S&P 500 performance), which give an indication about the client’s 
initial utility function 17(0, •). Then, the fund manager’s problem is to find a pair (17, vr*) with the 
given 17(0, •). Therefore, the question of finding large classes of forward performance processes U 
and the corresponding optimal portfolios vr* is of great importance. 


Assuming the filtration of the investor to be generated by a Brownian motion and her forward 
performance process U (t, x) to be an Ito process in t and twice continuously differentiable in x, one 
can show (see [MZ5] and [NT] for more details) that 17 is a solution of the fully nonlinear stochastic 
partial differential equation (SPDE) 


( 1 . 1 ) 


dU{t,x) = - 


1 \\Ux{t,x) X{t) + cr{t) a ^{t) (t,x)\\ 


17 XX (l) ^) 


■ dt -|- a{t, x)^ d>V(t). 


Here W = {W,B) is a standard Brownian motion that generates the filtration of the investor; 
VF is a Brownian motion to whose filtration the asset prices are adapted; cr is the corresponding 
volatility matrix of the asset prices and a~^ is its Moore-Penrose pseudoinverse; A is the market 
price of risk; a = (a'^, o'®) is a suitable stochastic process adapted to the filtration of the investor; 
and the superscript T denotes transpose. 


The forward SPDE (|1.1|) provides the analogue of the Hamilton-Jacobi-Bellman (HJB) equation 
that is associated with the classical optimization problems of expected utility from terminal wealth. 
As in the traditional setting, it is fully nonlinear and possibly degenerate. There are, however, 
fundamental differences between (ll.l|i and its classical counterpart. Firstly, (jl.ip is posed forward 
in time, which makes the problem ill-posed. Secondly, the forward volatility process a(t, x) is up 
to the investor to choose, in contrast to the classical case, where it is the mere outcome of the 
ltd decomposition of the value function process. The specification of the correct class of forward 
volatility processes is a very challenging problem, which remains open. 


So far three classes of forward performance processes have been exhibited in the literature: 1) 
time-monotone forward performance processes, that is: forward performance processes which are of 
finite variation in the time variable (see |MZ4j for more details); 2) homothetic forward performance 
processes, that is: forward performance processes whose dependence on the investor’s wealth x is 
of power form (see |NZ1| and [NTj for more details); 3) forward performance processes of factor 
form in complete markets (see [NT] for more details). These three types of forward performance 
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processes result from significant simplifications of the SPDE (11.111 in certain special cases: time- 
monotone forward performance processes are obtained by setting a = 0 and solving the resulting 
partial differential equation (PDE); homothetic forward performance processes result from choosing 
U to be of product form with power function dependence on x] and forward performance processes 
of factor form in complete markets are derived using a reduction of the SPDE (11.11) to a Hamilton- 
Jacobi-Bellman (HJB) equation that can be linearized in the complete market framework using the 
Fenchel-Legendre transform. 

As in |MZ5j and |NTj we consider forward performance processes of factor form (see Subsection 
II.II for the exact details), so that (jl.ip can be reduced to an (ill-posed) HJB equation. However, we 
consider the incomplete market case in which the HJB equation cannot be linearized by a simple 
transformation. Nonetheless, we are able to find explicit formulas for the leading order and first- 
order correction terms of the solution to such an HJB equation. These yield the leading order and 
first-order correction terms of the corresponding pair (U,tt*). We complement these results by an 
appropriate estimate of the remainder term. 

A similar expansion for the classical Merton problem in an incomplete market was given in 
[FSZj . In contrast to the Merton problem setup in |FSZ| . we face the additional difficulty of the 
HJB equation being ill-posed. In addition, no general estimates of the remainder were given in 
[FSZj . so that our approach provides new insights in the Merton problem setting as well. 

The following subsection describes our framework. 


1.1. Setting. We consider n tradeable securities whose prices follow the stochastic differential 
equations 

(1.2) dSiit) = Si{t)^ii{Y^{t),Y^{t))dt + Si{t)a^{Y^{t),Y^{t)fdW{t), i = 1, 2, ..., n 

and where Y^, Y^ are two observable real-valued stochastic factors. The factors are modelled by 
one-dimensional diffusion processes 

dY^{t) =5biY^{t))dt + V6K{Y^{t))dBi{t), 

(1.3) dW(t) = - ^{Y^it)) dt + ^ afY^it)) dB2{t). 

e ^/e 

We think of 6, e as being small positive numbers, so that Y^ should be thought of as a slow factor 
(e.g. a macroeconomic indicator) and Y^ as a fast factor (e.g. a fast mean-reverting stochastic 
volatility). Hereby, the noise W = (W, Bi, B 2 ) is jointly a {d + 2)-dimensional Brownian motion, 
VP is a d-dimensional standard Brownian motion, Hi, B 2 are one-dimensional standard Brownian 
motions, and the covariance structure is given by 

p^-t := {Wj,Bi){t), j = 1,2, d, 

p]t:= {Wj,B2){t), j = l,2,...,d, 

p^'U:= {Bi,B2){t). 


Since we allow for non-perfect correlation between the asset price processes and the stochastic 
factors, the market is in general incomplete. 


In our setting the forward performance SPDE (jl.ll) reads 

1 \\U,{t,x)X{Y^{t),Y^{t))+a{Y^{t),Y^{t))a{Y^{t),Y^{t))-\^{t,x)f 

(1.4) 2 Uxx{t,x) 

+ a{t, xY d>V(t). 
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Here (T(y^(t), y^(i)) = ,... is the volatility matrix of the stock 

price processes o, and 

X{Y^{t),Y^{t)) = {a{Y^{t),Y^it)f)-^f^{Y^{t),Y^{t)) 

is the market price of risk. The superscripts T and —1 denote transpose and Moore-Penrose 
pseudoinverse as before; and W is the standard Brownian motion obtained from the Brownian 
motion W by left-multiplication with a suitable constant matrix. 


In this paper we focus on solutions of (II.4p of factor form, namely on processes represented as 
(1.5) U{t,x) = V{t,x,Y^{t),Y^{t)), 

for some deterministic function V = V{t,x,yi,y 2 )- As discussed in the introduction of [NT] , such 
solutions are particularly natural from the economic point of view. Indeed, thinking of the forward 
performance process U as encoding the preferences of the investor on a set of trading strategies 
given the state of the world she observes and assuming that there are only finitely many quantities 
the investor keeps track of, it is natural to assume that the state enters her preferences through 
the corresponding finite number of factor processes. 

Assuming V £ applying Ito’s formula to V(t, x, Y^{t),Y‘^{t)), equating first the resulting 

martingale part with the martingale part on the right-hand side of (II.4p . and then the two bounded 
variation parts, one concludes that the function V{t,x,yi,y 2 ) is a classical solution of the HJB 
equation 

1 II Vx A -|- Cr a~^ (Yxy-i Vs Kp^ + Vxvi 'V P^) 11^ 

.6) Vt+A^/V-'-Yl - ^ =0. 

Vxx 

Here generator of the diffusion process (Y^,Y^). 

We also note that the initial condition U{0,-) for the SPDE (11.41) translates into an initial 
condition V{0,x,yi,y2) for the HJB equation (II.6p . so that the latter is posed in the “wrong” time 
direction and, in particular, one does not expect solutions to exist for all initial conditions or to 
depend continuously on them. In general, this ill-posedness is the main mathematical difficulty in 
dealing with forward performance processes. 


The main results of the paper (Proposition 14.11 Theorem 14.41) identify explicitly the leading 
order and first order correction terms of the solution V of (|1.6I) in the limit regime J 0, e 0. 
This allows to identify the leading order and first order correction terms of the corresponding pair 
(U,TT*) explicitly as well (see Proposition [531) • results are obtained under the following two 

assumptions. 


Assumption 1.1. (i) The range of left-multiplication by the matrix a is all of so that a a~^ 

is the d X d identity matrix. In particular, this implies n > d (and so the incompleteness in 
the model stems from the imperfectly correlated factors). Further, A is smooth in (yi,y 2 )- 

(ii) The initial condition does not depend on the factors, and we will write it as 14(0, x). Further 
F(0,x) is increasing and strictly concave in x. 

(iii) The function (14)(“^)(0, e“^) admits the representation 

(K)(-')(0,e-")= / ^!l^i.o(d^) + Co, 

Jr ^ 

for some non-negative finite Borel measure uq on M and constant Cq £ M, where (l/r)^”^^^)') 
is the inverse of the function 14 ( 0 , •). 


The latter condition is related to the ill-posedness of the initial value problem for the HJB 
equation (11.61) . and will turn out to be necessary for the leading order term of V to be well-defined. 
A class of possible initial conditions is given by F(0, x) = cix'^^^ > 0, C 2 £ (0,1). 
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Assumption 1.2. The process is positive recurrent with a unique invariant distribution /r. Clearly, 
the latter does not depend on the value of e (since a change in e corresponds to a multiplication of 
the generator of by a constant). 

1.2. Outline. To ensure that the main ideas are not obscured by cumbersome notation we first 
consider the cases where only the slow factor Y^ is present (“slow factor case”, Section [2]) or only 
the fast factor Y’^ is present (“fast factor case”, Section [3|). 

In the slow factor case we provide explicit formulas for the leading order and first order correction 
terms of V in Propositions 12.11 and 12.21 and justify the approximation of V by such in Theorem 
12.61 The corresponding results in the fast factor case can be found in Propositions 13.ll and 13.21 and 
Theorem 13.61 In Section 0] we consider the general case and give explicit formulas for the leading 
order and first order correction terms of V in Proposition 14.11 The corresponding remainder 
estimate can be found in Theorem 14.41 Finally, in Section [5] we give explicit formulas for the 
portfolios associated with our approximation (Definition 15.ip and explain in which sense they are 
approximately optimal fProposition l5.3p . 


2. Forward investment problem with a slow factor 


The first situation we consider is the slow factor case, that is when Hi and cjj in (11.2p depend 
only on Y^, and so V{t, x, yi,y 2 ) in (11.5p does not depend on y 2 . Moreover, to simplify the notation 
we write p for p® and y for yi throughout the present section. In view of Assumption ll.il and with 
these notations, the HJB equation (|1.6p becomes 


( 2 . 1 ) 


yt + 1^5<yYVyy+5b{y)Vy-- 


1 Wv^cKv) Y^xy'/bi^{y) p\? 


V,, 


= 0 . 


Here we aim to find an expansion of V of the form 

( 2 . 2 ) V = VY) + + 0(S) 

in the limit regime 5 j, 0. To this end, we will first derive expressions for and informally 
and then justify the resulting expansion in Theorem 12.61 below. 

2.1. Asymptotic Analysis. To obtain the leading order term we set 5 = 0 in m- 


(0)n2 


(2.3) 


,40) 1 l|A(y)f (K' O 

yi T7 


f(°) 

y XX 


= 0 . 


In addition, we endow the latter equation with the initial condition VYYo,x,y) = F(0, x). The 
resulting problem corresponds to taking the volatility coefficient in the forward performance SPDE 
(11.41) to be zero. This is precisely the case of time-monotone forward performance processes studied 
in |MZ4] . The formula for the solution F^^^ of (j2.3p can be therefore recovered directly from |MZ4[ 
Theorems 4 and 8]. 

Proposition 2.1 (Leading order term, slow factor). The solution F*-^) of the HJB equation (12.3p 
with the initial condition vY){f)^x,y) = F(0,x) admits the following representation in terms of uq 
and Cq in Assumvtion \l.l\f iii): 

(2.4) V^y{t,x,y) = u(^\\\{y)ft,x 
where u is given by 

(2.5) u{t,x) = — / e~^^ ^^^’Y~^ 2 hx[s,hy^\s,x)'jds + V{0,x), 

2 Jo 


h{t, x) = 


/ 

4R 




- 1 


uo{dz) + Co, 









6 


MYKHAYLO SHKOLNIKOV, RONNIE SIRCAR, AND THALEIA ZARIPHOPOULOU 


and denotes the inverse of h in the variable x. 

This follows from a transformation of (12.,‘ljl to the ill-posed heat equation and Widder’s Represen¬ 
tation Theorem of positive solutions to this equation f [Wil Theorem 8.1]). Both this transformation 
and Widder’s Theorem will be used to construct higher order terms of the expansion. The interpre¬ 
tation of (12.411 is that, at principal order, the forward performance measure is the complete market 
solution, but with the Sharpe ratio frozen to \{y). 

Next, we turn to the correction term in (12.211 . To obtain an equation for we plug 
R^*^^ -|- into (j2.ip and collect the terms of order To this end, we note the expansions 

II (k(») + (D)^ a + (f(»> + Vs VSKpf 


+ 


^(2RW Rji) ||Af + 2RW rW k + 0(6 ), 




r(0) 

^ XX 




Rii^ 

(Ri°V 


+ 0(5). 


The resulting equation for R^^^ reads 

^(1), 1 iiA(i/)f ^(1) 

. +2 ^yiO)y 


( 2 . 6 ) 


r(0) 

^ XX 


(0) T/0) T/(*^) 

-vi'^ = Vy)Hyfp 


r(0) 

^ XX 


and it is endowed with the initial condition R(^)(0,x,y) = 0, since the zeroth order term R^*^) has 
already satisfied the initial condition for R. 

Proposition 2.2 (Correction term, slow factor). The solution R^^^ of the PDE (I2.6jl endowed with 
the initial condition V^^'^{0,x,y) = 0 is given by 


(2.7) 


^ * s \\s \T yx^\t,x,y)v!^f{t,x,y) 

>{t,x,y) = -K{y)X{y) p -^-. 

2 Vfl\t,x,y) 


Proof. We start by introducing the change of variables 

(2.8) {t, y) := (^t, - log t, y'j , 

and set {t,^,y) = Vy\t,x,y), w^^\t,f^,y) = V^^\t,x,y). The latter functions are well-defined, 
since R(0) is strictly increasing and strictly concave in x, so that ^ is a strictly increasing function 
of X. 


By viewing the equation (12.3p as the “linear” equation 


(2.9) 


R 


(0) 


l|A(y)||^ 


R 


( 0 )' 


r(0) 
y XX 


Ri°^-||A(y)||' 


R 


( 0 )' 


r(0) 
y XX 


Ri°) = 0 


with coefficients depending on R^*^\ computing the derivatives 


(2.10) 

R/°) = ) 

(2.11) 

Rjo) = u;f ) 

(2.12) 



■ w 


( 0 ) 


R 


( 0 ) 


r(0) 

Vxx z X 


- 2 


r(0) 
r XX 

m 


R 

T/( 0 ) \ 2 
V XX 


R 


(0) 


+ w 


(0) 


R 


(0)x 2 


R 


( 0 ) 


R 


(0) 


r(0) 
y XX 


X 
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and plugging (12.101) . (|2.11l) and (|2.12l) into (12.9p . we obtain 


(0) , l|A(^/)f V, 


( 0 )' 


wl M 


y(0) 


w 


(0) ( ^xx 






( 0 ) 


+ w. 


( 0 ) 




( 0 ) 


y(0) 

^ XX 


-||A(y)|p 




(0)' 


y(0) 

^yxx , 


W 


( 0 ) 


T/(0) 
^ XX 




(0) 


= 0 . 


This readily simplifies to 

(2.13) 




in the coordinates 

A similar computation for shows that (12.6p transforms into 
(2.14) = K X{y)'^ pt X{y)'^ X'{y) 

in the new variables. Hereby, to obtain the right-hand side of (j2.14p we have relied on the following 
two considerations: 

(a) Note from (j2.5l) that u is differentiable in t and that ut is differentiable in x, which, with the 
assumed smoothness of A, implies the differentiability of in y. Then, by differentiating 
(j2.13p in y and rearranging one obtains 


= -Kyf^'iy) 


w 


( 0 ) 
««■ 


Moreover, the latter equation, endowed with the initial condition (0, ^, ?/) = 0, has the 
unique solution rc® = —t X{y)'^ A'(y) (note that the uniqueness of the solution is a conse¬ 
quence of Widder’s Theorem). In the original coordinates this solution reads 


(2.15) = -tX{yfX'{y) 

(b) In addition, by using (j2.1ip . (j2.12p one finds 


( 0)3 2 




y(0) 

y XX 


V, 


( 0 ) 


14 


X Vxx 


( 0 ) • 


(0) 

= - 


^(0) 
y XX 


y(0) 

yxx 


V, 


( 0 ) 


X Vxx '' X 


( 0 ) 


A combination of (a) and (b) gives the right-hand side of (|2.14l) . At this point, one can check that 
the unique solution of ()2.14p . endowed with the initial condition w^^\0,^,y) = 0 is given by 




= ^-^i^iy)>^iyy p>^iyy >^'iy)wf^\- 

Hereby, the uniqueness part of the statement follows by applying again the Widder’s Theorem. To 
obtain the proposition, it remains to change the coordinates back to (t, x, y) and use that 


= -tX{yfX'{y) 


(Hi°y 

which can be obtained from (12.151) by a differentiation in x. 


V 


(0) 




(0) 


X Vxx 


□ 


Remark 2.3. This result can be considered as the forward performance analog of [FSZp Proposition 
3.3] for the (backwards in time) Merton problem, but here the transformation (12.8p is crucial to 
reduce to the ill-posed heat equation for which Widder’s Theorem can be applied. In the traditional 
utility maximization problem, the correction term can be computed directly using commutation of 
certain operators and Black’s (fast diffusion) equation for the risk tolerance function, which is the 
approach used in [FSZj . Moreover in [FSZ] , only the case of a single stock is considered, and here 
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the analysis is done with multiple assets to obtain the correct coefficients in the approximations, 
which could not be obtained from just the one asset case with simple parameter replacements. 

We see from (j2.7p that depends on the slow factor through the parameters n and A frozen 
to the values K{y) and A(y), as well as A'(y). It is easily computed in terms of derivatives of the 
complete market forward performance measure We now give an additional representation 

of the correction term which has a natural interpretation in terms of the original forward 
performance problem. 


Proposition 2.4 (Natural parametrization of correction term, slow factor). The solution of 
the PDE (12.61) endowed with the initial eondition V^^\0,x,y) = 0, written as w^^\t,f,,y) in the 
eoordinates {t,f,,y) defined in (12.81) . admits the representation 


(2.16) 


w^^Ht,C,y)= [ w^^'>’%t,f,y)ds 
Jo 


where eaeh is the solution of the initial value problem 

(2.17) 


W- 




= 0, t > s 


t -r ^ 

with initial condition 

(2.18) = s it(y) A(y)^ p ||A(y)f y). 

In particular, each can be represented as 

(2.19) w^^^’^{t,^,y) = [ p'*’2^(dz) 

Jr 

with being a signed finite Borel measure on M. 

In the original coordinates, the same representation reads 

( 2 . 20 ) V’'^\t,x,y) = [ V^^'^’^{t,x,y)ds 

Jo 

where each is the solution of the initial value problem 




wliL* + 




(Oh 2 


y(o) 

^ XX 


with initial condition 


n Vx^\s,x,y)Vfy’{s,x,y) 

{s,x,y) = K{y)\{y) p -^-. 


40), 


vM\s,x,y) 


Remark 2.5. We remark that each of the processes 1/dh® (or, equivalently, can be viewed 

as an “auxiliary” forward performance process. These should be interpreted as the first order 
corrections that the investor makes at any given time s in reaction to the market conditions she 
observes. Furthermore, when making a projection of her future preferences from a time t onwards, 
the investor corrects her leading order forward performance criterion 17^*^) (or, equivalently, 
by aggregating all her previous first order corrections P^)’®, s G [0,t] (or, equivalently, tc^^^’®). 
We refer to as auxiliary to stress that it is neither a complete market forward performance 

measure, nor a solution of the full incomplete market problem. 
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Proof of Proposition We first recall that is a classical solution of the forward heat equation 
(I2.13p . Hence, the same is true for At this point, an application of Widder’s Theorem shows 
that the initial value problem (|2.17p . (j2.18p has a solution that exists for all t > s. In other words, 
each function s > 0 is well-defined. Since the forward heat equation with source (12.141) 

has a unique classical solution starting from the zero initial condition by Widder’s Theorem, the 
representation (I2.16P will follow once we establish that the right-hand side of (I2.16P is a classical 
solution of (I2.14p . This is the result of the following computation: 

= tK{y) X{y)'^ p\\X{y)f wf^\t,^,y) - (t,f,,y) ds 

= tK{y)X{y)'^p\\X{y)fwf^\t,^,y)- y) • 

Finally, (I2.19h follows applying again Widder’s Theorem [Wip Theorem 8.1], and the representation 
(|2.20p is the result of writing (I2.16jl in the original coordinates {t,x,y). □ 


The next theorem shows that, under appropriate assumptions, the error in the approximation of 
the true value function V by is indeed of order 6, as one would expect. To this end, 

we define the non-linear functional 


( 2 . 21 ) 


(0) T^(0) 




xy 


^(0) 

y XX 


+ 2 "^ 


1/2 
|2 xy 


2 ^ 


Vyy hVy , 


recall the change of coordinates {t,x,y) 1 —)• (t,f,,y) of (|2.8p . and set 


( 2 . 22 ) 


= r]\t,x,y). 


The bound on the approximation error can be then stated as follows. 


Theorem 2.6 (Remainder estimate, slow factor). Suppose that there exist do > 0 and T < 00 such 
that for all 6 G (0, do) the HJB equation (12.11) has a solution V G C'^’^’^([0, T) x (0, 00 ) X M) which 
is increasing and strictly concave in the second argument. Then, 

(i) the quantity 


(2.23) 


/ 

Jr 


E 

fc =0 


(- 1 ) 


k ^2k 


{2k)\2^ 


2k 


/ d \ P r P 

tk\^) Jq 2||a|G. dxdsdz. 


with as in (12.221) . is well-defined and finite for all d G (0, do), and 
(a) for every {t,x,y) G [0,T) x (0,oo) x M for which the limit superior 
(2.24) 

(_l)fc^2fc 


lim 

< 54.0 


/ 

Jr 


A 

2||A|t^( 


E 

fc =0 


(2A:)!2'= ||A||2fc t^ 


/ d \ 2^ ri r _ 

V^j Jo 2||X||^S dxdsdz 


is finite, the error bound 


lim d ^ 
(54-0 


V(t, X, y) - (t, x,y) - Vs {t, x, y) 


< 00 


(2.25) 


applies. If the limit superior in ()2.24l) is bounded above uniformly on a subset of [0, T) x 
(0, 00) X M, then the limit superior in (I2.25P is bounded above uniformly on the same subset 
of [0, T) X (0, 00) X M. 
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Remark 2.7. The meaning of condition (12.241) can be understood as follows. As explained in the 
proof of Theorem 12.61 below, the nonlinearity arising in the expansion of the HJB equation (12.ip 
is fed into an initial value problem for a backward heat equation through a source term in the 
new coordinates {t, y). The latter problem is severely ill-posed, with its solution operator rapidly 

magnifying the inverse Laplace modes (i.e. the analogues of Fourier modes for the inverse Laplace 
transform operator) of the source term. Therefore, in order to control the solutions uniformly for 
all small positive 5, one needs an a priori estimate on the inverse Laplace transform of the source 
term that is uniform for all small positive 6. The latter is precisely the content of condition (j2.24l) . 
In fact, the proof of Theorem 12.61 reveals that condition (j2.24l) is sharp in the sense of (j2.24p being 
equivalent to (|2.25p . 


Proof of Theorem 12.61 We start by expressing the HJB equation (j2.ip as 

r2 

11 A 11 “ 1 / “ ,— ^ »/ „ 1 / n ., 1 / 

(2.26) Vt - 


2 V^V^y 5 2 

p^ = -^ 

^XX ^XX ^ 


2 Fry 0 2 y _XhY 

T/ 9 ^ ^yy OOVy. 

* XX ^ 


Next, we write V = + VSQ, insert the latter expression into the left-hand side of 

(j2.26p . and expand in ^/6 using the elementary identity 


1 


1 


= i-^ ^ V , 

a + \/j b CL + Vs ab 


a< 0 , 6 <-^, 


Recalling that the functions and were constructed in such a way that all terms in (I2.26P 
on the orders of 1 and VS cancel out, and collecting the remaining terms we obtain after a lengthy 
but straightforward computation that 


(2.27) 




( 0 ) 


VSQt-V5\\XV^Q^ + Vs 


( 0)^2 


y(0) 

^ XX 


jvr) 

(Hi°^) 


Qa 


Sr]^. 


Next, we let Q := VS Q, recall the change of coordinates 

(L f., y) ■= (l - log 14 ® (L X, y) - t, y^ 

of (|2.8p . and define qitViV) '■= Q{t^cc,y). By the same computation as in the proof of Proposition 
[za the partial differential equation (|2.27p can be rewritten as 

(2.28) qt + ^ = Sff 

where v) = cc, y) as before. Then, Duhamel’s principle for the backward equation (12.281) 

implies that 

/TT^iiMi f VtV-X,y) e"2||A||"t dx = S [ [ V{sV-X,y)e~^^^^ dxds. 

V 27 rt||A|| Jr Jo v 27 rs||A|| Jr 

It follows that the right-hand side of the latter equation is in the domain of the inverse Weierstrass 
transform in the sense of |Wi2[ equations (5), (6)]. This, in turn, implies that the quantity in (I2.23p 
is well-defined and finite for all 5 £ (0,(5o)- Moreover, applying the inverse Weierstrass transform 
to both sides of (I2.29p . we see that, up to a multiplicative constant, the function ^ 1 —)• q{t,f,,y) 

is given by the quantity inside the absolute value in (j2.24p . The statement of the theorem is now 
immediate. □ 
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2.2. Example: Forward Performance Process of Power Type. We illustrate the results with 
the family of power utility forward performance processes. For a constant risk aversion coefficient 
7 G (0, oo)\{l}, we impose the initial condition for the HJB equation (II. 6 p to be 

P(0,x)=7'^f^. 

1-7 

This corresponds to Examples 16 and 18 in |MZ4| . 

We focus on an example where there is an exact solution to the problem with one factor volatility. 
Namely, we take a market model 

dS^{t) = Si{t) iii{Y\t)) dt + Si{t) ai{Y\t))'^dW{t), i = 1, 2, ..., n 
dY^{t) = 6im -Y^{t))dt + V5p ^Y^{t) dBi{t), 
where cri{y) and yi{y) are chosen such that 

A(y) = A ^/y, for some A G 

and the Brownian motions Wj are correlated with the Brownian motion Bi according to 
{Wj, Bi)(t) = pjt. We further suppose that ||/ 9 |p / This class of model was used, for 

instance, for a classical inhnite horizon consumption problem in |CV| . 


We insert the ansatz 

V{t,x,y) = j - g{t,y) 

1-7 

into the HJB equation (j2.ip to end up with 

2 

9 t + i^^ 0 ^y 9 yy + - y) + r ^/Jy( 5 lsJp) gy + ^ {\\H^y9 + ^ 0 ^y \\pf = 0, 

where F = and the initial condition is g{0,y) = 1. Next, we make the transformation g = 
with the choice q = which after a short computation leads to the linear equation 

+ \si3^y ^yy + {6{m - y) + F v^/3 A^p) '^y + ^ II Af y = 0 

with initial condition y) = 1. Looking for a solution of the form exp(Ai(t)y + A 2 {t)) we hnd 
that Ai, A 2 need to solve the ODEs 


a; = ^5p^Al + (VjF/ 3 A^p - 6) A, + ^ ||Af, 

A!^ = 5m Ai 

with zero initial conditions. Let a± be the roots of the quadratic corresponding to the right-hand 
side of the Riccati equation for Ap. 

^ + (VJF/3 A^p - J)a + = 0 . 

2 2q 

We assume that the model parameters are such that both roots are real and a_|_ > 0. This is the case 
when 7 G (0,1), A'^p < 0, and F q (A^p)^ > ||A|p; or 7 > 1, A'^p > 0, q < 0, and F q (A^p)^ > ||A|p; 
or 7 > 1 and q > 0. In any of these cases, we then obtain the solution 

1 — / 2 
Ai{t) = a_ ^ ^ and A 2 {t) = 5m log 
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where A is the square root of the discriminant of the quadratic above. We, in turn, deduce that 
since a+ > 0 , the solution 

(2.30) V{t,x,y) = 7 ^'^— exp {q {Ai{t) y + A 2 (t))) 

is well-defined for all f > 0 . 


We can now check condition (j2.24l) of Theorem l2.6l to show the convergence of our approximation 


1/^°^ {t, x,y) + V6 (t, X, y) = 7 '^ 




1-7 




2r 7 
-1 - 

7 1-7 




where the expression for follows from Proposition 12.11 with vq being a Dirac mass at 7 ^ and 
Co = 7 , and the expression for is a consequence of Proposition 12.21 


In (|2.21l) . the term in the parentheses in the hrst summand turns out to be zero, because both 
V and are multiples of ''', and the term in parentheses in the second summand is of order 
Vs, because can be obtained from V by setting <5 = 0, and V is smooth in S. Therefore 
r/ is a product of x^~'^ and a function of (t, y) depending smoothly on 6 (which follows from the 
smoothness of Ai, 242 ), and hence, V ™ (12.22j) is a product of and another function of {t,y) 
depending smoothly on 6. An explicit computation of the inverse Weierstrass transform in (j2.24jl 
reduces the expression there to 


lim 

54,0 


f g^r (t s)r 2 iiji F(^s,y,5) ds 

Jo 


for some function F depending smoothly on 6. In particular, the latter limit superior is hnite and 
so the conditions of Theorem 12.61 are satisfied in this example. 


3. Forward investment problem with a fast factor 


The second situation we consider is the fast factor case, that is when /ij and in (II.2|) depend 
only on and so V{t, x, yi,y 2 ) in (|I.1>P does not depend on yi. To simplify the notation, we write 
p for and y for y 2 throughout this section. With these notations and in view of Assumption lI.il 
the HJB equation (II. 6 |) becomes 


(3.1) 


P + 


®(y)^ i{y) 1 11^ ®(y) pIP 

^ ^ - 2 


= 0 . 


Our goal here is to hud an explicit expansion of the solution V to (13.ip of the form 


(3.2) V = + Ve + 0(e) 

in the limit regime e 0 . As in the previous section, we will hrst derive formulas for and 
informally, and then justify the resulting expansion by means of a suitable remainder estimate. 

To hnd we plug (13.2p into (13.11) and collect the leading order terms (namely, those on the 
order of e“^) to get 

(3.3) hhlv-m +7fe)V7(0) _ ^„(„)2||^||2h^ 

Yxx 

Note that we can satisfy (j3.3p by choosing as a function of t and x only. As we explain below, 
the exact choice of P^®^ will be pinned down by the lower order terms in the expansion of (13.ip . 
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To proceed, we plug (13.21) into m, and collect the terms of order e We obtain 

^XX ^xx \^xx ) 

Choosing to be independent of y as we noted earlier, the latter equation simplihes to 

Clearly, the above equation can be satisfied by choosing to be a function of (f, x) only. As 
with , the exact choice of will result from considering lower order terms in the expansion 
of (l3dH . 

Next, we insert the extended expansion 

V = + e + O(e^) 

into (13.11) in order to find the terms of order 1. This results in the equation 

,3.4) 


y(0) 

y XX 

Hereby, we have used the fact that and do not depend on y. The next proposition shows 
that there is a unique choice of such that (0, x,y) = V (0, x) and the equation (13.4p . viewed 
as an elliptic partial differential equation for has a solution. 

Proposition 3.1 (Leading order term, fast factor). There is a unique function such that 
V^^^0,x,y) = 17(0, x) and the equation (13.4p possesses a solution . With 


A = 


l|A(i/)rMdy) 


\ 1/2 


sueh a function admits the representation 

(3.5) V^^\t,x,y) = V^^\t,x) = u{X^t,x) 

where u is given by 


1 

u{t,x) = —- J /i*^“^)(s,x)^ ds + 17(0,x), 

f _ 1 

h{t,x)= / - no{dz) + Co. 

Jr z 

Here h^~^'> is the inverse of h in the variable x, and uq o.iT'd Cq were introduced in Assumption 
\l.lY iii). 

Proof. We start by integrating (13.4p with respect to the invariant distribution y. of Assumption [T21 
Since 17 does not depend on y and 

^ + i{y) ) /i(dy) = 0 

(due to the invariance of y), we obtain 


(3.6) 


(0) A2 (Hi°y 


- 

* 2 


1/® 
^ XX 


= 0 . 


We easily conclude using |MZ4[ Theorems 4 and 8]. 


□ 
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From ()3.5I) we observe that gives the complete market forward performance measure with 
constant Sharpe ratio A, where the appropriate averaging has been identified by the asymptotic 
analysis. 

To obtain we will expand the HJB equation (|3.1I) up to order which however requires 
further information on As a first step, we subtract from the equation ()3.4I) its averaged version 
(13.61) to get 

|2 _ p (Vi0))2 


(3.7) 


aivY 


yi?+7{y)Vy^^^ = 


We introduce the notation 
(3.8) 


r 

<P{y) = / 

Jo 




||A(yH«))f-A^ 


y(0) 

y XX 


ds, 


where denotes the fast factor process, solution of the SDE (II.3p . but with e = 1. Then (see, for 
example, [FPSSl Section 3.2, p. 94]), the solution of (|3.7p admits the stochastic representation 

(3.9) + C{t,x), 

where C{t,x) is a function that does not depend on y. 

We can now expand (|3.ip up to order to obtain 

+ Kfi^-^A(y)^(^y«A( 2 /) + (/)'(y)(^^^^) a{y) 


+ = 0 . 


Averaging this equation with respect to the invariant distribution p of Assumption 11.21 we obtain 

further 

(3.10) 


r,"> + 


Tr 


( 0 )' 


y(0) 
, y XX 


Ty(0) 

y(i) _ ^2 - 

''xx ( 0 ) '^X 

^xx 


V, 


(0) 


y XX 


(yi°y 

vM' 


Y'iy) a(y) Hyf yi^y )) p = o. 


which is the desired partial differential equation for Since satishes the initial condition 

for V, we endow p.lOl) with the initial condition x, y) = 0. 

Proposition 3.2 (Correction term, fast factor). The unique classieal solution of the partial differ¬ 
ential equation (I3.10p with the initial eondition V^^\0,x,y) = 0 is given by 


T^(0) / 

(3.11) V^^Yt,x,y) = V^^Yt,x) = t-^ 
with (j) as in dMl). 

Proof. We introduce a new space variable 

(3.12) 


(yi°y 


(Y{y) a{y) X{yf p{dy) p, 


A2 


^=-iogyio)-yt. 

Since is strictly increasing and strictly concave in x for any given t, f, is a strictly increasing 
function of x and we may define 


w 


( 0 ) 


{t,f,):=V^^\t,x) and w^^\t,f,):=V^^\t,x). 
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Following the calculations in the proof of Proposition 12.21 gives 

(3.13) inf ^ i P u;® = 0. 

A similar computation starting from (|3.10l) shows that the transformed correction term satisfies 
the forward heat equation 

(3.14) inf^ ^ ^ K^y)^ p > 

with the initial condition in^^fO,^) = 0. At this point, it is easy to check (using (I3.13|) i that 

(3.15) in^^) = t wf^ 4>'{y) a{y) A(y)^ p{dy)^ p 

satisfies (I3.14p with the desired initial condition. Changing back to the original coordinates we 
easily obtain (13.lip . 

The uniqueness part of the proposition follows from the uniqueness of the solution of the Cauchy 
problem for the forward heat equation (see |Wi[ Theorem 8.1]). □ 

Remark 3.3. This result can be considered as the forward performance analog of [FSZp Proposition 
2.7] for the (backwards in time) Merton problem, but here the transformation (I3.12p is crucial to 
reduce to the ill-posed heat equation for which Widder’s Theorem can be applied. The differences 
highlighted in Remark 12.31 apply here too. 

Next, we give an additional representation for the correction term which has a natural 
interpretation in terms of the original portfolio optimization problem. 

Proposition 3.4 (Natural parametrization of correction term, fast factor). Let be the 

functions from Provositions \3.1[ \3.S\ written in the coordinates 

■= (i,-logF® - 

Then: 

(i) admits the representation 

(3.16) = f ^) ds 

Jo 

where each is the solution of the initial value problem 

(3.17) -k y = 0, t > s, 
with initial condition 

(3.18) 0 = wf^ {s, 0 (p'iy) a{y) X{y)'^ hidy)^ p. 

It can be therefore represented as 

(3.19) wdh^t,0= [ u^^\dz), 

Jm. 

with being a suitable signed finite Borel measure on M. 
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(ii) In the original coordinates, the same representation reads 

yW(t,x)= [ V^^^’^{t,x)ds 
Jo 


where each is the solution of the initial value problem 


^ ^ ^ = 0, t> 


( 0 ) 






y(0) 

^ XX 


with initial condition 




V^^\s,x) 

Vx^{s,x) V Vxx\s,x) ) X 


(j)'(y) aiv) X{yf fi{dy)^ 


P- 


Remark 3.5. The quantities (or, equivalently, of Proposition [33] should be interpreted 

in the same way as “auxiliary” forward performance processes as their analogues in the slow factor 
case. We refer to Remark 12.51 above for more details, but point out that the asymptotic analy¬ 
sis identifies the constant vector [J^4>'{y) a{y) X{y)'^ ij{dy)) p as the principal correcting effect of 
stochastic Sharpe ratio, in the case of nonzero correlations p. 


Proof of Proposition \3.Ji\ We first recall from the proof of Proposition 13.21 that is a classical 
solution of the forward heat equation (|3.13l) . Hence, is a solution of the same equation and, 
therefore, the solutions s > 0 of (I3.17p . (I3.18P are well-defined and given by 


= wf^{t,C) (y) a{y) X{y)'^ p{dy)^ p. 


Therefore the right-hand side of the representation (I3.16P is equal to the right-hand side of (13.151) . 
and, thus, (I3.16P immediately follows. Moreover, the representation (13.191) is a direct consequence 
of Widder’s Theorem. Finally, part (ii) of the proposition can be either established in the same 
way as part (i), or by changing to the original coordinates {t,x) in (|3.16l) . (|3.17l) and (I3.18p . □ 


We conclude this section with the appropriate remainder estimate. Specifically, we will show 
that the error in the approximation of the true value function V by is of the order 
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e. To this end, we introduce the non-linear functional 
(3.20) 

+ ^ l|A(y)f ( 4i|) 7^ 

+ Kyf vP (A(y) + a{y)pV^l'>) (y - 

- i l|Afe)lP ri») (K - r<") - f.W) (- ' ) 


1 T/l'^) 1 1 

+ ^ A(?/)^ Q;(y) p (V,y - e Fj^)) + — — || A(y) V, - A(y) -7 e"^/^ a(y) p | 

1 \\v.Ky) + v^y^My)pf \ f^A^) ||A(^)f 

\ \ 2 / \ 2 


-7 


^ 3/2 


(1) , l|A(y)||- 


, / 


y(l) _ 


Vn 


( 0 ) 


14 


( 0 ) 


"l \(yffvPA(y) + ^'(y) 


a(y)p 


Here is defined through (13.91) . and we note that the value of does not depend on the 
choice of the constant C{t,x) in (13.91) . We also set := (t, — log t, y), and let 

fl%t,i,y) = p'^{t,x,y). 


Theorem 3.6 (Remainder estimate, fast factor). Suppose that there exist cq > 0 and T < oo such 
that for all e G (0,eo) the HJB equation (13.ip has a solution V G C'^’^’^([0, T) x (0,oo) x M) which 
is increasing and strictly concave in the second argument. Then, 

(i) the quantity 



(_l)fc^2fc 

{2k)\2^ t^ 



X,y)s 2 . dxdsdz 


(defined via (I3.20p and the paragraph following it) is well-defined and finite for all e G (0, cq), 
and 

(a) for every {t,x,y) G [0,T) x (0,oo) x M for which the limit superior 


(3.21) 


lim 

etO 


r 2 ^ 

e-^Y. 

k=0 


k ^2k 


i-irz 


d A 


{2k)\2^t^ \Ai) 


2k 


/7^'< 


s,^-X,y)s 2. dxdsdz 


is finite, the error bound 


(3.22) 


lim 

40 


-1 


V(t, X, y) - (t, x)-yfe V^^'> {t, x) 


< oo 


applies. If the limit superior (|3.2ip is bounded above uniformly on a subset o/[0, r)x(0, oo)xM, 
then the convergence in (j3.22p is uniform on the same subset of [0, T) x (0, oo) x M. 


Remark 3.7. Condition (I3.21j) is of the same form as condition (I2.24h . and the detailed interpretation 
of the latter given in Remark 12.71 applies here as well. 


Proof of Theorem \3.(A We proceed as in the proof of Theorem l2.6l Specifically, we insert the ansatz 
V = 1/^°^ R® -l-e^ Q ™to the HJB equation (13.ip and expand the resulting 
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equation in the powers of The terms V^‘^\ and were chosen in such a way that 

all terms on the orders of 1 and cancel out. At this point, a tedious but straightforward 

computation relying on the elementary identity 


1 = 1 _ fl/2 ^ 

a + e ^/2 6 a a^ + e^/^aft’ 


a < 0, b < —e a 


allows to compute the terms of order e and leads to 


( 0)12 


(3.23) 


^Qt + 


{viS'f 


Qxx ^ 


I|A(^)|PF. 


( 0 ) 


■ Qx = erf 


where Q := e~^ {V — +eQ and 77 '^ is dehned according to (13.201) . 

One can then conclude the argument by repeating the steps in the proof of Theorem 12.61 making 
the change of coordinates 


(T y) ■■= (t, - log 



in (13.231) . and combining Duhamel’s principle for the resulting equation with the formula for the 
inverse Weierstrass transform given in [Wi2| . □ 


Remark 3.8. Consider the example of Section lT^ with the explicit solution reparametrized according 
to (5 = e~^, and take the corresponding approximation for a fast volatility factor as considered in 
this section. Then, by direct computation (/)(y) = ||A|p(?/ — m) and, therefore, our approximation 
is 

+ Ve||A||2/3mA^p,r2 7T' 

1 — 7 1 ~ 7 

The corresponding remainder r]‘' has a complicated dependence on e. However, an explicit Taylor 
expansion in e of the functions Ai, A 2 in the formula for V of ()2.30l) (with 6 replaced by e~^) shows 
that 7 *^ is given by a product of x^~'^ and an order one in e function of t, y. Consequently, fj'^ is a 
product of and a function of {t,y) that is of order 1 in e. The argument of Section YI?2\ yields 
that the conditions of Theorem 13.61 are satisfied in this example and its fast factor approximation. 


4. Multiscale forward investment problem 

We combine our approaches to the forward investment problems with slow and fast factors 
to analyze the multiscale forward investment problem described in Section 11.11 We consider an 
expansion for V{t,x,yi,y 2 ) in equation (ll.6jl of the form 

V = + Vs + Ve + 0(S + e) 

in the limit regime <5 0, e 0. We first give the general results and, in turn, explicit formulas for 

the case of power utilities in Section 14.21 

4.1. First Order Approximations. It is convenient to define: 

(4-1) A(?/i) = ||A(?/i,y 2 )||^/^(d 2 / 2 )^ 

(4.2) Cl, 0 ( 1 / 1 ) = (p®)^ A(yi, ^ 2 ) p(di/ 2 )^ Vvi), 

Co,1(2/1) = (p'^)^ (^ J ^ X { yi , y 2 ) 4 > 1/2(2/112/2) “(2/2) p(d2/2)^ , 


(4.3) 

















FORWARD PERFORMANCE PROCESSES 


19 


where 

(4.4) 


/■OO r 9 - 1 

4>{yi,y2)= E ||A(yi,y^(s))|| - A^(yi) | y^(0) = ^2 ds, 
Jo 


and denotes the fast factor process, solution of the SDE (11.31) . but with e = 1. 

The following proposition gives explicit formulas for the leading order term and the Hrst 
order correction terms ydh) y(o,i)^ 

Proposition 4.1 (Explicit formulas, general case). (i) The leading order term admits the 
representation 

V^°\t,x,yi,y 2 ) = V^^\t,x,yi) = u {}?{yi)t,x) 

where u is given by 

u{t,x) = — / e~^^ x)) ds + 1A(0, x), 

2 Jo 


h{t,x) = [ 

4r 


^zx-\z'^t _ Y 


vo{<Jz). + Cq. 


Here h^ is the inverse of h in the variable x, and pq and Co were introduced in Assumption 
\l.lV iii). 

(a) The slow scale correction term jg ly 


(4.5) 


Vxx 


and admits the natural parametrization 

(4.6) V^^’^\t,x,yi)= /VW’^’^(t,x,yi)ds, 

Jo 

where each is a solution of the initial value problem 


il), 5 ,s ^ / vP \ ] yWAs = 0, t>s, 


( 0 )' 




with initial condition 


y(0) 
, ^XX 




( 0 ) 


X,») = CM 

vP{s,x,yi) 


(Hi) The fast scale correction term ig giyen by 


(4.7) 


V^'^’^Ht,x,yi) =tCo,i{yi) 




( 0 )' 


y(0) 

y XX 


(u'°y 
h? /1 


The function 1/(^4) admits the natural parametrization 

V<^°’^\t,x,yi) = /VW’^’^(t,x,2/i)ds, 
Jo 

where each is a solution of the initial value problem 






+ 


>?{yi 


V, 


( 0 )' 


y(0) 

y XX 


- \\yi) 


V, 


( 0 )- 


y(0) 

y XX 


e4),^,* = 0, t>s, 
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with initial condition 


= Co, 1 ( 2 / 1 ) 


Vx°\s,x,yi) 

V^x\s,x,yi) 


/ yi°^(g,x,2/i)^ 

V vM\s,x,yi) 


X 


Remark 4.2. The quantities of Proposition 14.11 should be interpreted in the same 

way as their analogues in the single factor cases. We refer to Remark l2.5l for more details. However, 
we highlight that the analysis identifies the following reduced parameters: A( 2 / 1 ) in (14.Ih . the Sharpe 
ratio root mean square-averaged with respect to the fast factor and frozen at the value 2/1 of the 
slow factor; Ci, 0 ( 2 / 1 ) in (14.2p . which has the effect of the correlation p® between the slow factor and 
equity returns; and Cq, 1 ( 2 / 1 ) in (j4.3p . which has the effect of the correlation between the fast 
factor and equity returns. 


Proof of Proposition \4-l\ We start with the proofs of parts (i) and (hi). To this end, we insert 
(5 = 0 into the HJB equation (11.61) , employ Assumption 11.11 and then proceed as in the proofs of 
Propositions EH ESI and Eai The arguments from therein can be repeated directly by replacing 
A(2/) by A( 2 /i, 2 / 2 )- In particular, and determined via an averaging of the equations 


(4.8) 
and 

(4.9) 


h»> - 


2 (Vjo))'^ 


y(0) 

y XX 


+ Cy^ ^ 0 




( 0 , 1 ) 


+ 




( 0 )' 


y XX 


a + ap\+ Cy, p(0’3) = 0 

A V \^xx Z X / 


with respect to p(d 2 / 2 )- Here e ^Cy 2 is the generator of the fast factor that is 

C ?/2 dy 2 y 2 +7(2/2)51/2) 

and the terms PPi, are the ones appearing in the expansion of the solution to 

(jl.6l) with (5 = 0. In particular, subtracting from (14.81) its averaged version, we obtain the expression 


1 

p(°’^^(t,x, 2 /i, 2 / 2 ) = -^<(>( 2 / 1 ) 2 / 2 ) ^ "^( 0 ) + C(t,x, 2 /i) 
where C{t,x,yi) is a function that does not depend on 2 / 2 , and cf was defined in (14.4p . 


It remains to prove part (ii). To this end, we again employ Assumption 11.11 and insert the ansatz 
into the HJB equation (jl.6p . Collecting the terms of order y/5 in the resulting 

equation we get 


(4.10) 

^(1) 


P® (/)'«P 




y(0) 

y XX 


+ 




^( 0 ) 

^xy2 

y(0) 
r XX 


+ ■ 


Aw 4" + ^ 


a V ' ' 
y2V2 


ye 


We can now write pi^i = -|- y^p(^d) _)_ £^( 1 , 2 ) expand equation (|4.10p in powers of e as 

in the proof of Proposition 13.21 By doing so, we conclude that and 1/(^4) ggg chosen as 

functions independent of 2 / 2 - 
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Moreover, can be determined by averaging the equation 


(4.11) + 


1 




( 0 )' 


y(0) 

y XX 


y(b0) _ 


4". 


(0) 


T/(0) T/{0) 

y xv^ yx 


^(0) 
y XX 


v-M _ „ AV + A, rM = 0 


y(0) 

y XX 


with respect to /i(dy 2 )- The averaged equation is endowed with the initial condition 
2 ;, _ Q pgg solved explicitly by means of a transformation to an ill-posed 

backward heat equation as in the proofs of Propositions 12.21 and 12.41 This gives the explicit formula 
for and its natural parametrization. □ 


Remark 4.3. Formulas (|4.6p and (|4.7p can be considered as the forward performance analogs of 
those appearing in [FSZp Section 4.1] for the (backwards in time) Merton problem. The differences 
highlighted in Remark 12.31 apply here too. 


We now complement the explicit formulas for the leading order terms p(0d) gy 

convergence theorem justifying the approximation of the true value function V by the function 
+ \/eP^^’^^. We will need the non-linear functional whose lengthy formula we 
give in AppendixWe also set (t, yi, 2 / 2 ) := (T — log t, yi, 1 / 2 ), and let 

(4.12) fj^’^{t,C,yi,y2) := y'^’"(t,x,yi,y2). 


Theorem 4.4 (Remainder estimate, general case). Suppose that there are 6 q > 0, cq > 0, and 
T < 00 such that, for all (5, e) £ (0,5o) x (OjCo)? the HJB equation (II.6p has a solution V £ 
P' 1 . 2 , 2,2 T) X (0, 00 ) X R^) which is increasing and strictly concave in the second argument. Then, 

(i) the quantity 


{S + e)-^ 



- (_l)fc^2fc 

^ (2fc)!2*^ t^ 



- X,yi,y2) S 2. dydsd^. 


with defined in (j4.12p . is well-defined and finite for all {d,e) £ (0,5o) x (0, eo), and 
(ii) for every {t,x,yi,y 2 ) £ [0,T) x (0,oo) x R^ for which the limit superior 


(4.13) 

lim (5 -|- e)“^ 
54 . 0,40 



(_l)fc^2fc 

{2k)\2^ t^ 



X,yuy2)s zsdydsdz 


is finite, the error bound 


(4.14) lim (5-|-e) ^ 

54.0,e4.0 


V{t,x,yi,y2) -V^°\t,x,yi) -VdV^^’^'>{t,x,yi) - y/eV^^'^\t,x,yi) 


< 00 


applies. If the limit superior (I4.13p is bounded above uniformly on a subset of [0, T) x (0, 00 ) x 
R^, then the convergence in (I4.14p is uniform on the same subset o/[0,T) x (0, 00 ) x R^. 


Remark 4.5. Condition (14.131) is of the same form as condition (|2.24l) . and the detailed interpretation 
of the latter given in Remark 12.71 applies here as well. 


Proof of Theorem \4.4\ We proceed as in the proof of Theorem 12.61 More specifically, we plug 
p = p(o) _|_ p(i.O) 4- y^p(Od) _|_ Q HJB equation (II.6p and apply a Taylor’s expansion 

to the resulting equation in y/5 and y^. Hereby, we use the elementary identity 

1 _ 1 b 

a + b a + ab 
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and the definitions of and to eliminate the ter m s of 

orders 1, and ^/e. The remaining equation then reads 

(«5) + 

\ Vxx ) \ ^xx ) 

where rf'’" is defined prior to the statement of the theorem. One can now conclude by repeating 
the steps in the proof of Theorem 12.61 namely by making the change of coordinates 

(4.16) {t,C,yi,y2) ■= (t,-logV^^'> -^^^t,yi,y2'j 

in ()4.15p . and combining Duhamel’s principle for the resulting equation with the formula for the 
inverse Weierstrass transform given in [Wi2| . □ 

4.2. Power utility example. We illustrate the results with the family of power utility forward 
performance processes. For a constant risk aversion coefficient 7 £ (0,oo)\{l}, we impose the 
initial condition of the HJB equation (II. 6 j) to be 

F(0,x)=7^f^. 

1-7 

Similar to Section Ea we have the following explicit solution for the constant parameter value 
function V^^\t,x,yi) in part (i) of Proposition 14.11 

(4.17) V^'^\t,x,yi) = u (X‘^(yi)t,x) where u(t,x) = ^'^— -and T = 

1-7 7 

which can be verified by taking the measure z/q to be a Dirac delta centered at 7 “^ and the constant 

Co = 7. 

From (14. 5 p in Proposition 14.11 we compute 

= ^t^Cifi{yi)T^X{yi)\'{yi)V^^\t,x,yi), 

and from (14.71) we obtain 

Then, the three-term approximation to the forward performance value function is given by 

(4.18) 

y{t,x,yi,y2) = (^1 + ^^t‘^Ci^o{yi)T‘^X{yi)X'{yi) + y/etCo,iiyi)T^'^ V^^\t,x,yi) + 0{d -Fe), 

where is given explicitly in (14.171) . 


5. Approximately optimal portfolio 

In this last section we define the portfolio associated with our approximation and establish its 
approximate optimality. 


Definition 5.1. Let the feedback portfolio function vr^ he given by 


7r~ = -{a ) 


(5.1) 


\ VXX 


T/(0) T/(liO) _ T/(lh) 

V XX V X Vx r XX 


(0)n2 




+ Vi 


T/(0) T/(0,1) _ T/(0) Ta(0,1) 
^ XX ^ X * X y XX 

M\2 




/T I T\-l s y^y! , T\-l f 4‘y2 f {y^ \ 

-y'Snia ) p ^ + ^e{a ) pa^A—^\ 

y XX y XX y y XX 


X 
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where and are as in Proposition \4-l\ and cj) is given in (sai). The formula dSU) is 

obtained by recalling that the nonlinearity in the HJB equation (HSI) results from the optimization 
problem 

sup (^{XVx + VdKp^ Vxyi + -^a Vxy 2 f {<y T^) + ^ Vxx (fT7r)^(o-^)^, 

replacing V by its expansion ^+ y/d + y^e y(o>i) ihg formula for the corresponding opti¬ 
mizer n*, and applying Taylor’s expansion in 6, e to the result. We refer to tt~ as the approximately 
optimal portfolio which is justified by the next proposition. 

Remark 5.2. One can use formula (I4.18P for the value function in the case of the power forward 
performance to compute the approximately optimal portfolio. We omit the lengthy expression here. 


Proposition 5.3. Suppose that there are 6 q > 0, cq > 0, and T < oo such that, for all 
(5, e) € (0,(5o) X (0, eo), the HJB equation ()1.6I) has a solution G T) x (0, oo) x R^) 

which is increasing and strictly concave in the second argument. Then, the value process 
V^’^{t, {t),Y^{t),Y^{t)) satisfies a SDE of the form 

(5.2) Vd K dBi (t) + Vyf a dB 2 {t) + dW{t) + 0^’" dt 

with drift coefficient 0*^’^ = Q^’^{t, (t), Y^{t),Y‘^{t)). If, in addition, the limits superior resulting 

from (|4.13p by replacing fj^’^ with any of pfy^ finite, then 

(5.3) lim ((i + e)“^ |0'^’''(t,x,?/i,y2)| < oo. 

o4..0,e4..0 

In other words, the performance of the portfolio 7r~ fails to fulfill the martingale criterion of opti¬ 
mality only by a bounded variation term of order 5 + e. 


Proof. The SDE (15.2p is obtained by combining 

dX^(t) = p{Y\t), W(t))^7r~(t) dt + a{Y\t),Y^{t))^Tr-{t) dW{t), 

the definition of 7r~ (Definition [5T]), and Ito’s formula. Writing -\-y/dV^^'^^ -\-^/eV^^’^^ + 

and recalling the HJB equation (II.6p we see that the drift coefficient 0^’'" is a linear combination 

of (5x^ Q% 1 , Qxy 2 , and Qxx, with the coefficients being uniformly bounded in d, e wherever the 

limits superior lim Vx’^, lim 14™, I, lim hm-W are finite. 

54.0,ejo (54.0,eto ' ' (54-0,e4.0 ' ' (54-0,40 V44 


Next, recall the change of variables (I4.16p . write q^''’{t,f,,yi,y 2 ) for (^'^’'^(t,x,?/i, 2 / 2 ), and recall 
from the proof of Theorem 14.41 that q^’^ satisfies the backward equation 

IIXI 12 

(5.4) 


S,e I 
Qt + 


S,e ~S,e 

= V’ 


with zero initial condition. Differentiation of this equation in Duhamel’s principle for the resulting 
equation, and the formula for the inverse Weierstrass transform given in |Wi2j yield that q^’'" (and 
therefore also Qx^) admits a control of the form (15.31) at a point if condition (I4.13P holds for 
at that point. Repeated differentiation in f gives the same result for q^’^ (and therefore also Qxx) 
and as well as for q^'^^ and 


Lastly, differentiate the equation (15.411 in yi to get 
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By linearity — t solves 


Qt + 


■% = 



Hence, admits a control of the form (15.31) wherever condition (14.131) is satisfied by , whereas 
q^’^^ has already been controlled before. The desired bound on Qxy^ then follows from the bounds 
on qfy-^, qf^, and q^’^. To finish the proof it remains to estimate Qxy 2 in a similar manner and to 
combine the estimates on Qx'^, Qxyi, Qxy 2 , and Qxx- □ 

Remark 5.4. In the case that only a slow factor is present (that is, in the setting of Section [2|) a 
statement analogous to that of Proposition 15.31 holds for the portfolio function 


TT 


= -{a^)-^X 


V, 


(0) 


y(0) 

y XX 


+ ^5 


yiO) Ty(l) _ T/(0) T/(l) 
y XX y X y X y xx 

(0)\2 




^ - V6K{a'^) ^ p 




where and are as in Propositions 12.ll and 12.21 respectively. 

We easily deduce that a similar result holds for the case that only a fast factor is present (that 
is, in the setting of Section [3]) for the portfolio function 


7r~ = — (cr^) ^ A 




( 0 ) 


y(0) 

yxx 


+ Ve 


T/(0) Ty(l) _ y(0) T/(l) 

yXX yX yx yxx 


\ . r, T,-l 4 >y 

'' y XX ^ y XX 


(viff 

where and P*-^) are as in Propositions I3.ll and 13.21 respectively, and (/>(y) is given in (|3.8I) . 


6. Conclusion 

We have provided a convergent approximation for forward performance processes in a multifactor 
incomplete markets model, as well as for the corresponding optimal portfolio. Our approach is 
based on a perturbation analysis of the corresponding ill-posed HJB equation. The principal term 
in the approximation results from the appropriately averaged problem, whose solution is known 
from Widder’s Theorem. The correction terms for fast and slow volatility factors can be computed 
explicitly in terms of this leading order term. 

We have also given explicit calculations in the case of power utility. The ease of the formulas 
provided in more general cases, as well as conditions for convergence, should allow future work 
to develop the financial implications of forward performance processes in realistic market environ¬ 
ments. 
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Appendix A. Expression for in Section H] 


: = 




(Er))^ + 


(^J0,l))2 + ^ _ yiO) _ ^^yil,0) _ ^vjO,l)y 


+ ^ {V.y, - - V~6V^lf^ - ^ + <5 kA'^ y Ejo) 

+ K A'^y Ejo) IIAf Eji-O) ^jo,i) ^ ^ ^ ^(o) yim 


Vs 


+ < 5^/2 ^ X^pS yilfi) y(l,0) + 5 ^ ^ yil,0) y(0,l) ^ ^ ^ _^T^f ^jRO) y^^^ 

+ VreK yW + 6 Vet, A^y e a A^y 

+ a A^/ E,,, + 5^/2 ^2 ||^s||2 yio) y(i,o) ^sVV ||y f E® Ejo^^) 

+ ^ a K Ei°) + 5^/2 ^ ^2 ||^S||2 yjl,0) yjO.l) + ^ Q, ^ Vxy 2 

ye ve 

^VSaK^fV (^x,i - eW - ^/5Eil;0) - ^/iE4V0 

^ {Vs \ VV yVS\ Ei^’O) + 5 K y E® + <53/2 Ei'f) + 5 Vi K y E^V^ + ^ ) 

+ A (E, - E^ - V5 Ei^’O) - V 

^ {Vs \ Eji’O) + Vi A EjOA) + V5 ic V E W + <5 A V Ej^f) + V^ k V Ej^f) + ^ « / K,.) 
+ V5 It A^V (E, - VV - Vs Ei^’O) - V {Vxy, - EiS - V^yf^ - V V.?4) 

2 VjeW eWa + ||a ||2 viEjo) vV^ + V^it aV VfS 


VIrx Eis^ 


+ 1 a Ejo) E,,, + II A||2 E® (E, - E^) - VJ - V vV^^) 
ye 


vV^Y (E,,-EiV) 


(0)n2 


-A: 


y _ Ym ^v -e E(2) _ V/2 E® _ V e^/2 y(i-2)^ _ 


^ \VxxY ^ 

Here Ai V^‘^\ At and At E*^^’2) are defined through (j4.8p . (I4.9p and (14.lip , respectively. 
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